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Abstract 

^^ I The projectile fragmentation reactions using ^^Ni & ^'^Ni beams at 140 MeV/n on targets ^Be 

& ^^^Ta are studied using the canonical thermodynamical model coupled with an evaporation 









code. The isoscaling property of the fragments produced is studied using both the primary and the 
secondary fragments and it is observed that the secondary fragments also respect isoscaling though 
the isoscaling parameters a and /3 changes. The temperature needed to reproduce experimental 
data with the secondary fragments is less than that needed with the primary ones. The canonical 
model coupled with the evaporation code successfully explains the experimental data for isoscaling 
for the projectile fragmentation reactions. 
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I. INTRODUCTION 

Projectile fragmentation reaction is used extensively to study the reaction mechanisms in 
heavy ion collisions at intermediate and high energies. This is also an important technique 
for the production of rare isotope beams and is used by many radioactive ion beam facilities 
around the world. The fragment cross sections of projectile fragmentation reactions using 
primary beams of "^^Ca, ^^Ca, ^^Ni and ^"^Ni at 140 MeV/nucleon on ^Be and ^^^Ta targets 



have been measured at the National Superconducting Cyclotron 



aboratory at Michigan 
al has been used to cal- 



State University [ll]. The canonical thermodynamical model(CTM) 
culate some of these fragment cross sections [sI]. In the present work, an evaporation code has 
been developed and has been coupled with the canonical thermodynamical model. CTM 
coupled with this secondary decay code is then used to analyze the isoscaling data from 
the projectile fragmentation reactions. The lighter fragments produced from these reactions 
exhibit the linear isoscaling j^-|6| phenomena and our model calculation also strongly sup- 
ports this observation. The secondary fragments (produced after applying the evaporation 
code on the canonical model) also exhibit isoscaling like the primary fragments from the 
fragmentation reaction [TJ but the temperature needed to reproduce the experimental data 
with the secondary fragments is lower than that required by the calculation with the pri- 
mary fragments. The isoscaling parameters a and /3 as obtained in the present work from 
the model calculations agree closely with those obtained from the experimental data. These 
parameters as obtained from the secondary fragments are lower in magnitude than those 
obtained from the primary ones. This effect is also seen in the dynamical models |8| though 
the reduction is much more there. The isoscaling behaviour displayed by the fragments 
produced in the projectile fragmentation reactions and the effect of evaporation on it in the 
framework of the HIPSE model has been discussed recently in [9|. 

This paper is structured as follows. First we describe the canonical model briefly in 
Sec. II. In the same section we also present the main features of the evaporation code which 
is used to calculate the secondary fragments. In Sec. Ill we present the results. The effect 
of sequential decay on the distribution of the isotopic fragments is discussed. The isoscaling 
phenomena as displayed by the primary as well as the secondary fragments is also described 
in this section and are compared with the experimental data. In Sec. IV we present the 
summary. 



II. THE STATISTICAL MODEL 

In models of statistical disassembly of a nuclear system formed by the collision of two 
heavy ions at intermediate energy one assumes that because of multiple nucleon-nucleon 
collisions a statistical equilibrium is reached. Consequently, the temperature rises. The 
system expands from normal density and composites are formed on the way to disassembly. 
As the system reaches between three to six times the normal volume, the interactions between 
composites become unimportant (except for the long range Coulomb interaction) and one 
can do a statistical equilibrium calculation to obtain the yields of composites at a volume 
called the freeze-out volume. The partitioning into available channels can be solved in the 
canonical ensemble where the number of particles in the nuclear system is finite (as it would 
be in experiments). In the next subsection we describe the canonical model. 

A. The canonical t her mo dynamical model 

In this section we describe briefly the canonical thermodynamical model. Assume that the 
system with Aq nucleons and Zq protons at temperature T, has expanded to a higher than 
normal volume and the partitioning into different composites can be calculated according 
to the rules of equilibrium statistical mechanics. In a canonical model, the partitioning is 
done such that all partitions have the correct Aq, Zq (equivalently Nq,Zq). Details of the 



B; 



implementation of the canonical model can be found elsewhere [2| ; here we give the essentials 
necessary to follow the present work. 

The canonical partition function is given by 

"^■^ 

Here the sum is over all possible channels of break-up (the number of such channels is 
enormous) which satisfy Nq = J^^ ^ ^i,J and Zq = ^J ^ ^jj; ujj is the partition function 
of one composite with neutron number I and proton number J respectively and rij^j is the 
number of this composite in the given channel. The one-body partition function uj^j is a 
product of two parts: one arising from the translational motion of the composite and another 
from the intrinsic partition function of the composite: 

ui,j = ^{2nmTf'A'/' X ziA^nt) (2) 



Here A = I + J is the mass number of the composite and Vf is the volume available for 
translational motion; Vf will be less than V, the volume to which the system has expanded 
at break up. We use Vf = V — Vq , where Vq is the normal volume of nucleus with Zq protons 
and A^o neutrons. In this calculation we have used a fairly typical value V = GVq. 

The probability of a given channel P{njj) = P{nQi,niQ,nii rijj ) is given by 

P{n,,) = ^U^ (3) 

The average number of composites with / neutrons and J protons is seen easily from the 
above equation to be 

{ni,j) = ^i,J—7^ (4) 

The constraints Nq = Y,I >^nj j and Zq = Y, J ^^i j can be used to obtain different looking 



but equivalent recursion relations for partition functions 10| ■ For example 



Qno,Zo - tt 51 ^^i,jQno-i,Zo-j (5) 

These recursion relations allow one to calculate Qno,Zo 

We list now the properties of the composites used in this work. The proton and the 
neutron are fundamental building blocks thus zifi{int) = zo^i{int) = 2 where 2 takes care 
of the spin degeneracy. For deuteron, triton, ^He and ^He we use zi^j{int) = (2s/,j + 
1) exp{—f3Eij{gr)) where /3 = l/T,Eij{gr) is the ground state energy of the composite 
and {2si^j + 1) is the experimental spin degeneracy of the ground state. Excited states for 
these very low mass nuclei are not included. For mass number A = 5 and greater we use 
the liquid-drop formula. For nuclei in isolation, this reads {A = I + J) 

The derivation of this equation is given in several places [2|, [ll| so we will not repeat the 
arguments here. The expression includes the volume energy, the temperature dependent 
surface energy, the Coulomb energy and the symmetry energy. The term - — represents 
contribution from excited states since the composites are at a non-zero temperature. 

We also have to state which nuclei are included in computing Qno,Zo (eq-(17)). For J, J, 
(the neutron and the proton number) we include a ridge along the line of stability. The 



liquid-drop formula above also gives neutron and proton drip lines and the results shown 
here include all nuclei within the boundaries. 

The long range Coulomb interaction between different composites can be included in an 



approximation ca 
scheme set up in 



led the Wigner-Seitz approximation. We incorporate this following the 



ll|. 



B. The evaporation code 

The statistical multifragmentation model described above calculates the properties of the 
collision averaged system that can be approximated by an equilibrium ensemble. Ideally, 
one would like to measure the properties of excited primary fragments after emission in 
order to extract information about the collisions and compare directly with the equilibrium 
predictions of the model. However, the time scale of a nuclear reaction(10~^°s) is much 
shorter than the time scale for particle detection (10~^s). Before reaching the detectors, 
most fragments decay to stable isotopes in their ground states. Thus before any model 
simulations can be compared to experimental data, it is indispensable to have a model 
that simulates sequential decays. A Monte Carlo technique is employed to follow all decay 
chains until the resulting products are unable to undergo further decay. For the purposes of 
the sequential decay calculations the excited primary fragments generated by the statistical 
model calculations are taken as the compound nucleus input to the evaporation code. Hence, 
every primary fragment is decayed as a separate event. 

We consider the deexcitation of a primary fragment of mass A, charge Z and temperature 
T. The succseesive particle emission from the hot primary fragments is assumed to be the 
basic deexcitation mechanism. For each event of the primary breakup simulation, the entire 
chain of evaporation and secondary breakup events is Monte Carlo simulated. The standard 
Weisskopf evaporation scheme is used to take into account evaporation of nucleons, d, t, 
He^ and a. The decays of particle stable excited states via gamma rays were also taken into 
account for the sequential decay process and for the calculation of the final ground state 
yields. We have also considered fission as a deexcitation channel though for the nuclei of 
mass < 100 its role will be quite insignificant. The process of light particle emission from 
a compound nucleus is governed by the emission width F,^ at which aparticle of type u is 



emitted. According to Weisskopf's conventional evaporation theory 



12-], the partial decay 



width for emission of a light particle of type v is given by 
gmo^ {E* -Eo- K) 



r, = ^L^'- '- ^ exp(2 Ja«(E* - E^ - K) - 2JapE*) (7) 

TT^/l ttR 

Here m is the mass of the emitted particle, g is its spin degeneracy. Eq is the particle sep- 
aration energy which is calculated from the binding energies of the parent nucleus, daughter 
nucleus and the binding energy of the emitted particle and the liquid drop model is used to 
calculate the binding energies. The subscript z/ refers to the emitted particle, P refers to the 
parent nuclei and R refers to the residual(daughter) nuclei, ap & a^ are the level density 
parameters of the parent and residual nucleus respectively. The level density parameter is 
given by a = A/16MeV~^ and it connects the excitation energy E* and temperature T 
through the following relations. 



E* = apTl 
{E* -Eo- K) = aRTl (8) 

where Tp & Tr are the temperatures of the emitting(parent) and the final (residual) 
nucleus respectively. V^ is the Coulomb barrier which is zero for neutral particles and non- 
zero for charged particles. In order to calculate the Coulomb barrier for charged particles of 
mass A> 2 we use a touching sphere approximation J13l| . 



K = uf" " ^'^^'^ 1/3 for A. > 2 

n{A'J' + {Ap - A^f} 

iZp - l)e2 
= -p^ — tor protons (9j 

nAj,"" 

where r, is taken as 1.44m. 

(To is the geometrical crosssection (inverse cross section) associated with the formation of 
the compound nucleus (parent) from the emitted particle and the daughter nucleus and is 
given by (Tq = ^ri?^ where. 



R = ro{(Ap - A,f'^ + A,i/3} for A, > 2 



where tq = 1.2 fm. 

For the emission of giant dipole 7-quanta we take the formula given by Lynn|l4| 



Pp{E 
with 



^i = -T^f depn{E*-e)f{e) (11) 



.. . 4 1 + ^ e^ NpZp Tgs' 

^^^^ SirmnC^hc Ap {TGey + {e^ - E^y ^ ' 

with K = 0.75, and Eq and Tq are the position and width of the giant dipole resonance. 
For the fission width we have used the simplified formula of Bohr- Wheeler given by 

r/ = ^exp(-5;/Tp) (13) 



where Bf is the fission barrier of the compound nucleus given by 15 1 



Bf{MeV) = -lAOZp + 0.22{Ap - Zp) + 101.5. (14) 

Once the emission widths are known, it is required to establish the emission algorithm 
which decides whether a particle is being emitted from the compound nucleus. This is 



done [l6| by first calculating the ratio x = r/rtot where Ttot = fi-/^tot, ^tot = J2u^u and 
u = n,p, d, t, He^, a, 7 or fission and then performing Monte-Carlo sampling from a uniformly 
distributed set of random numbers. In the case that a particle is emitted, the type of the 
emitted particle is next decided by a Monte Carlo selection with the weights T^/Ttot (partial 
widths). The energy of the emitted particle is then obtained by another Monte Carlo 
sampling of its energy spectrum. The energy, mass and charge of the nucleus is adjusted 
after each emission. This procedure is followed for each of the primary fragment produced 
at a fixed temperature and then repeated over a large ensemble and the observables are 
calculated from the ensemble averages. . The number and type of particles emitted and the 
final decay product in each event is registered and are taken into account properly keeping 
in mind the overall charge and baryon number conservation. 

III. RESULTS 

First we will show our calculations for ^^Ni on ^Be reaction and ^'^Ni on ^Be reaction. In 
the model, the target imparts a certain amount of energy to the projectile transforming it to 



a projectile like fragment (PLF) with a temperature. This excited PLF will then expand and 
form composites during the expansion. The partioning of the PLF into different composites 
is done by the rules of equilibrium statistical mechanics in a freeze-out volume. We consider 
production of different isotopes from the statistical breakup of the dissociating system. If 
< riij > is the average number(multiplicity) of composites with i neutrons and j protons, 
then the cross-section for this composite is <j{i-,j) = C < Uij >, where C is a constant not 
calculable from the thermodynamic model. It depends upon the dynamics that are outside 
the scope of this model. To be able to compute < Uij > we need to know the mass and 
charge of the PLF and its temperature. The source sizes adopted for this calculation are zero 
order guesses. It could be sometimes smaller or greater depending on the diffusion from the 
target. For ^^Ni or ^^Ni on ^i?e which is a small target the choice of the mass and the charge 
of the PLF is limited. It can be slightly less than that of the projectile to as large as that 
of the projectile plus ^Be, the last being the case when the much larger projectile swallows 
the small target and drags it along retaining PLF features. Similarly we have some limits 
on energy imparted(this fixes the temperature). This energy can be small or upto the upper 
limit. The upper limit is given by the case of projectile swallowing Be and all the energy 
transforming into internal excitation(no part going into collective flow). In the canonical 
calculation, the dissociating system is taken to be ^^Ni +^ Bc^Nq = 35, Zq = 32) and for 
the other reaction the dissociating system is taken to be ^^Ni +^ Be{NQ = 41, Zq = 32). All 
composites between drip lines are included as detailed in Sec.IIA with the highest values of 
N, Z terminating at A'"o, Zq. The temperature is taken to be 5.8 MeV for both the reactions. 
Fig. 1 displays the isotopic distribution for Z= 12 (magnesium) produced from both the 
the reactions. The dashed lines correspond to the distributions of the primary fragments 
while the solid lines correspond to the distributions after sequential decay. As expected, 
the more neutron rich system with Nq/Zq = 1.28 (right panel) produces more neutron rich 
isotopes than the neutron deficient system with Nq/Zq = 1.09 (left panel). In all cases, the 
primary distributions are much wider and more neutron rich than the final distributions. 
The peak positions of the isotopic distributions of both the primary and the secondary 
fragments coincide in case of the neutron deficient system as seen from the left panel of 
the figure. In case of the neutron rich system(right panel) the peak of the distribution of 
the secondary fragments has shifted to the left with respect to that of the primary. The 
experimental isotopic distributions(solid squares with error bars) agree much more with the 
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final results obtained after secondary decay than with the primary distributions. The width 
and peak position of the isotopic distribution after the secondary decay agrees very well with 
the experimental data. The model also successfully reproduces the rapid fall in crosssection 
for large neutron number. 

We will now discuss the results about isoscaling. It is observed from the experimental 
data[17] that the light fragments emitted from the ^^Ni and ^^Ni systems exhibit the linear 
isoscaling behaviour represented by the equation 

R21 = Y2{N, Z)/Yi{N, Z) = Cexp{aN + (5Z). (15) 

where the isoscaling ratio R2i{N,Z) is factored into two fugacity terms a and /3, which 
contain the differences of the chemical potentials for neutrons and protons of the two reaction 
systems. Y2{N,Z) refers to the yield of fragment(N,Z) from system 2 which is usually 
taken to be the neutron-rich one and Yi{N,Z) refers to the same from system 1. C is a 
normalization factor of the isoscaling ratio. It is observed from our model that both the 
primary as well as the secondary fragments exhibit isoscaling. Fig. 2 shows the isoscaling 
results for Ni on Be system for the primary fragments. The ratio -R21 is plotted as function 
of the neutron number from Z =6 to Z= 13 in the left panel whereas the right panel displays 
the ratio as function of the proton number Z from N=8 to N=15. It is seen that the primary 
fragments exhibit very well the linear isoscaling behaviour for the lighter fragments over a 
wide range of isotopes and isotones. The lines in the figures are the best fits of the calculated 
R21 ratios(open triangles) to Eq.l5. They are essentially linear and parallel on the semi log 
plot . 

Fig. 3 displays the isoscaling results for the secondary fragments. The open triangles 
are the results obatined from our model while the solid squares with the error bars are the 
experimental ratios. We have shown the isoscaling results for the even Z and odd Z isotopes 
in two separate panels for the sake of clarity. While comparing with the results of the primary 
fragments in Fig 2, it is evident that the isoscaling is valid for a limited range of isotopes for 
the secondary fragments as compared to the primary ones. When the isotopes away from 
the valley of stable nuclei are considered, the trends for the secondary fragments are not as 
clearly consistent with the isoscaling law as are the trends of the primary distribution. One 
can conclude from this that isoscaling is approximately valid in the case of the secondary 
fragments. The lines in the figures are the best fits of the calculated R21 ratios which agree 



closely with the experimental data. The temperature required to reproduce isoscaling data 
with the primary fragments is about 8 MeV[7| whereas that required for the secondary 
fragments is 5.8 MeV. This decrease was already predicted in one of our earlier papers 7| 
in Sec. 9. It has also been found out by Ono et al.[8] that the effect of secondary decay is 
to decrease the isoscaling parameter a by about 50%. This is indeed what emerges from 
our calculations after including secondary decay code with the canonical model though the 
amount of reduction is less as compared to the dynamical model. Tha value of a as obtained 
from the fits of the primary fragments (left panel of Fig. 2) is 0.713 whereas that obtained 
from the secondary fragments is 0.580 which is much closer to the experimentally obtained 
value for q;[17| equal to 0.566. 

In Fig. 4 we have also plotted the isoscaling ratios for different neutron number values 
as function of the proton number Z and thereby calculated the other isoscaling parameter /3 
from them. The value of /3 as obatined from the linear fits of the primary fragments( right 
panel of Fig. 2) is -0.849. For the secondary fragments, the value of /3 from our model is 
-0.634 whereas the experimentally obtained value is -0.621. The value of /3 also decreases in 
the case of secondary from the primary ones as in the case of a. The results obtained after 
the sequential decay matches very well with the experimental ratios. As for the isotopes in 
Fig 3, for the isotones also it is seen that the linear isoscaling is valid for a limited range in 
case of secondary fragments as compared to the primary ones(right panel of Fig. 2). 

We now turn to the case where the target is ^^^Ta. One can consider the extreme limit 
which is target independence and by which we mean that A^o, -^o refers to simply to the 
case where just the projectile is the disintegrating system. The more likely scenario where 
some matter has diffused to or from the target has many possibilities. In principle, the 
target could shear away some material from the projectile leaving a PLF which is a fraction 
of the projectile. For a peripheral collision this is less likely than the alternative of the 
projectile picking up part of nuclear matter from the tail region of the much larger target. 
The amount of nuclear matter curved out of Ta will be small (for the disintegrating system to 
retain PLF characteristics) but other than that not much can be said and an integration over 
the different possibilities might be essential. Keeping this limitation in mind, we compared 
with different possible scenarios and the case with projectile plus 10 neutrons and 8 protons 
from the target, i.e, projectile +^^^0 yielded the best results. It is also seen that target 
independence gives worse results than this case. The agreement of our model calculations 
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with experimental data is pretty good in this case. The temperature used for this reaction 
is 6.2 MeV. 

Fig. 5 and Fig. 6 are similar to Fig. 3 & Fig. 4 except for the fact that they are for the 
case where the target is ^^^Ta instead of ^Be. The theoretical slope a as obtained from Fig. 
5 from the fit of the secondary fragments is 0.459 while the experimentally obtained slope 
is 0.432. The straight line fit to the calculated points matches nicely with the experimental 
ratios. The value of (3 as obtained from our model from the slopes in Fig. 6 is -0.489 whereas 
the experimentally obtained value is -0.487. Thus we find that the values of the isoscaling 
parametrs as obtained from the secondary fragments for both the targets agree quite well 
with the experimental values as can be seen from Table 1. 



Target material 


Isoscaling parameters 


primary 


secondary 


experiment 


^Be 


Q^ [Zmin = 6, Zmax = 13) 


0.713 


0.580 


0.566 


^Be 


l3{Nrmn = 8,N^ax = l5) 


-0.849 


-0.634 


-0.621 


^^^Ta 


Oi\Zmin = 0) ^max = 13) 


0.619 


0.459 


0.432 


^^^Ta 


P Nmin = 8, Nmax = 15) 


-0.682 


-0.489 


-0.487 



TABLE I: Best fit values of the isoscaling parametrs a and /3 for the two targets '^Be and ^^^Ta. The 
values obtained from the slope of the primary and secondary fragments as well as the experimental 
values are tabulated. In the second column the range of .Z^ or A^ values used to calculate the 
parameters are indicated. 



IV. SUMMARY 



This work deals with the developing of the sequential decay code and successfully coupling 
it with the canonical thermodynamical model in order to compare the properties of the 
secondary fragments with the experimental data. The width, peak position and rapid fall 
in cross-section of the isotopic distribution of the secondary fragments matches well with 
the experimental data. The main purpose is to examine the effects of sequential decay 
on the phenomenon of isoscaling. The secondary fragments also shows isoscsling and the 
the isoscaling parameters calculated from the secondary fragments matches closely with the 
experimental data. The temperature required to reproduce the experimental data with the 
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secondary fragments is less than needed by the primary ones. We finally conclude that 
the canonical thermodynamical model can explain the isoscaling properties of the lighter 
fragments produced in the projectile fragmentation reaction. 
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FIG. 1: Experimental cross sections of magnesium isotopes (squares with error bars) compared 
with theoretical results: primary fragments ( open stars joined by dotted line) and secondary frag- 
ments(open triangles joined by solid line). The left panel is for the reaction ^^Ni on ^Be while the 
right panel is for ^^Ni on ^Be reaction. The temperature is 5.8 MeV for both the reactions. 
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FIG. 2: Ratios(i?2i) of multiplicities of primary fragments of producing the nucleus (A^, Z) where 
reaction 1 is ^^Ni on ^Be and reaction 2 is ^'^Ni on ^Be. The left panel shows the ratios as function 
of neutron number N for fixed Z values from 6 to 13, while the right panel displays the ratios as 
function of proton number Z for fixed neutron numbers from A^ = 8 to 15. The lines drawn through 
the theoretical points(open triangles) are best fits of the calculated ratios. The temperature used 
for both the reacions is 5.8 MeV. 
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FIG. 3: Ratio of multiplicities of the secondary fragments of producing the nucleus {N, Z) where 
reaction 1 is ^^Ni on ^Be and reaction 2 is ^^Ni on ^Be compred with the ratios of the experimental 
cross sections of the same two reactions. The left panel shows the even Z isotopes while the right 
panel shows the results for the odd ones. The lines drawn through the theoretical points(open 
triangles) are best fits of the calculated ratios. The experimental points are shown by solid squares 
with error bars. The temperature used for both the reacions is 5.8 MeV. 
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FIG. 4: Same as Fig. 3 except that the ratios are plotted as function of the proton number Z for 
fixed neutron numbers. The left panel shows the results for the even neutron numbers while the 
right ones show those for the odd ones. 
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FIG. 5: Same as Fig. 3 except that here the target is ^^^Ta instead of ^Be. The temperature 
used for both the the reactions is 6.2 MeV. 
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FIG. 6: Same as Fig. 4 except that here the target is ^^^Ta instead of ^Be. The temperature 
used for both the reactions is 6.2 MeV. 
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